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Abstract: We present in vivo volumetric images of human retinal micro- 
circulation using Fourier-domain optical coherence tomography (Fd-OCT) 
with the phase-variance based motion contrast method. Currently fundus 
fluorescein angiography (FA) is the standard technique in clinical settings 
for visualizing blood circulation of the retina. High contrast imaging of 
retinal vasculature is achieved by injection of a fluorescein dye into the 
systemic circulation. We previously reported phase -variance optical 
coherence tomography (pvOCT) as an alternative and non-invasive 
technique to image human retinal capillaries. In contrast to FA, pvOCT 
allows not only noninvasive visualization of a two-dimensional retinal 
perfusion map but also volumetric morphology of retinal microvasculature 
with high sensitivity. In this paper we report high-speed acquisition at 125 
kHz A-scans with pvOCT to reduce motion artifacts and increase the 
scanning area when compared with previous reports. Two scanning schemes 
with different sampling densities and scanning areas are evaluated to find 
optimal parameters for high acquisition speed in vivo imaging. In order to 
evaluate this technique, we compare pvOCT capillary imaging at 3x3 mm 2 
and 1.5x1.5 mm 2 with fundus FA for a normal human subject. Additionally, 
a volumetric view of retinal capillaries and a stitched image acquired with 
ten 3x3 mm 2 pvOCT sub-volumes are presented. Visualization of retinal 
vasculature with pvOCT has potential for diagnosis of retinal vascular 
diseases. 
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1. Introduction 

At present, the most widely used technique for examining the circulation of the retina in 
clinical settings is fundus fluorescein angiography (FA). This is a procedure in which the 
fundus is photographed through a bandpass filter following excitation of an extrinsic 
fluorophore, fluorescein sodium (a compound that fluoresces after absorbing short- 
wavelength light) as it is injected into the blood stream through a vein. In an effort to replace 
this invasive procedure, one commercial system (Retinal Functional Imager, Optical Imaging) 
has been developed to image retinal vasculature without use of any extrinsic contrast agents 
using a fundus camera. In this paper we concentrate on optical coherence tomography (OCT) 
[1,2] methods that have the benefit of offering perfusion maps or flow with depth information 
in the human retina [3-8]. The phase-variance based motion contrast technique for Fourier- 
domain OCT (Fd-OCT) [9] has been implemented successfully in visualizing vasculature of 
the zebrafish [10], ocular circulation of the mouse [11], and human retinal vasculature [12]. 

Visualization of the retinal vasculature and microcirculation is essential for clinical 
diagnostics and treatment monitoring of retinal vascular diseases. The primary function of 
capillaries in the eye is to supply nutrients, exchange respiratory gases, and to excrete 
byproducts of metabolic reactions in the retina. Ocular vascular diseases transform the 
morphology of vessels and the formation of capillary networks, angiogenesis and vaso- 
occlusions caused by vascular endothelial cell aging, and damage associated with lifestyle- 
related illnesses such as diabetes mellitus, hypertension, and heart disease. In the human 
retina, there are generally three specific layers with capillary networks. The major arterial and 
venous branches are present in the nerve fiber layer or the ganglion cell layer. In addition, two 
capillary branches of the intermediate plexus and the deeper plexus are located in between the 
inner plexiform layer and the inner nuclear layer as well as between the outer plexiform layer 
and the outer nuclear layer, respectively [13]. In histology, the lumen diameter of human 
retinal vessels has approximately -100 |im of arteries and -200 |im of veins close to the optic 
disk and -15 |im of microcapillaries in the deeper retina [14]. In a clinical environment, 
fundus FA generates vasculature maps of the retina and blood stream flow for identification of 
arteries and veins as well as hemorrhage. Additionally, commercial system fundus FA 
supports relatively easy operation, fast capture, and high-resolution fundus images. In rare 
instances, however, FA may cause adverse reactions in patients who are hypertensive or 
elderly [15]. Moreover, patients may experience nausea, vomiting, or acute hypertension 
during or after fundus FA imaging. 

In order to measure blood flow and dynamics, optical Doppler tomography [16] was 
developed to measure velocity profiles of fluid flow by detection of Doppler shifts. Doppler 
OCT demonstrated in vivo flow imaging of the rodent skin [17] and the directional color map 
of biological tissue [18]. Since then, a number of in vivo vasculature imaging techniques have 
been reported using variations of OCT including: Doppler Fd-OCT with the histogram-based 
method [19], color Doppler Fd-OCT [20], Doppler Fd-OCT with the bulk motion removal 
method [21], optical coherence angiography [22], resonant Doppler flow imaging [23], joint 
spectral and time domain OCT [24], optical micro -angiography [25], speckle variance OCT 
[26], Doppler optical frequency domain imaging [27], and dual-beam Doppler OCT [28,29]. 
Phase-variance optical coherence tomography (pvOCT) identifies regions of motion between 
consecutive B-scans, which allows mapping of micro-vasculature without quantitative flow 
information [12]. 

In this paper we present contrast agent free images of in vivo human retinal vasculature 
using phase-variance based motion contrast OCT. High-speed image acquisition at 125,000 
A-scans/s, compared to previously reported 25,000 A-scans/s [12], allows a dramatic decrease 
of motion artifacts and increases the field of view that can be imaged within a few seconds in 
a single volume acquisition. To visualize capillaries in two- and three-dimensions, two 
scanning patterns are used with different sampling densities and scanning areas including 3x3 
mm 2 and 1.5x1.5 mm 2 . The processed data show depth-resolved retinal micro-vascular maps 
for different field sizes and total acquisition times. Furthermore, two-dimensional phase- 
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variance data over a selected depth region as well as three-dimensional representation of the 
vasculature are presented and compared with fundus FA images. 

2. Materials and methods 

2.1. Phase-variance optical coherence tomography instrument 




Fig. 1. Schematic of the Fourier-domain OCT instrument. The CMOS (Complementary Metal 
Oxide Semiconductor) camera is driven by a Lab VIEW operating program on a workstation 
computer through a frame grabber. Subjects can see a green dot from the fixation target 
through the dichroic beamsplitter. SLD: Superluminescent diode. 

Figure 1 is a schematic of the Fourier-domain OCT system with the CMOS 
(Complementary Metal Oxide Semiconductor) line detector [30] used to acquire data 
presented in this paper. A superluminescent diode (Superlum, Ireland) with an 855 nm central 
wavelength and a full width half maximum of 75 nm is used as a light source. The axial 
resolution is 4.5 Lim in the retina (refractive index, n = 1.38). The lateral resolution of the 
system at the retina is estimated at 10 Lim to 15 Lim depending on the subject's ocular 
aberrations. A standard fiber-based OCT instrument is operated in Michelson interferometer 
configuration, connecting the fiber coupler (OZ Optics, Ontario, Canada) with 90/10 splitting 
ratio. 10% of the input light is sent to the sample arm and the power of the imaging beam at 
the cornea is 350 llW. The remaining 90% of the light is sent toward the reference arm where 
it is back reflected by the mirror and passed to the detector arm where it interferes with the 
light back scattered from the sample. The reference arm length matches the sample arm length 
which is set by the position of the imaging sample. The resulting interference in the form of 
the spectral fringe pattern between the sample and the reference arm is detected in the custom- 
built spectrometer. The spectrometer consists of a 75 mm focal length collimating objective, a 
1200 1/mm diffraction grating, and a 150 mm objective which focuses the light onto the 
CMOS detector (spL4096- 140km, Basler). The maximum line rate of the detector is 140 kHz 
at 4096 pixels and 10x10 Lim (dual line) pixel size. Configuration setups of the CMOS camera 
are 4-tap, 12-bit, 80-MHz camera link clocks; the line averaging mode, a 125 kHz line rate 
and 2048 active detector pixels. Subject's gaze is stabilized using a custom-built fixation 
target reflected through a dichroic beamsplitter (FF662-FDi01, Semrock). The fixation target 
is composed of an 8x8 LED (light emitting diode) matrix which emits 516 nm (green) light 
and an electronic controller to switch fixation light position. During acquisition with 
Lab VIEW software (National Instruments), 3 to 5 seconds depending on scanning pattern 
settings, the raw data are stored in the workstation (HP xw8600, 3. 2 -GHz dual processors) 
with the rapid speed frame grabber (PCIe 1429, National Instruments). 

2.2. Image acquisition and data processing 

Fd-OCT data sets were acquired in the Vision Science and Advanced Retinal Imaging 
laboratory (VSRI) at the University of California Davis Medical Center on a 59 -year-old 
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healthy subject with normal ocular media. Written informed consent approved by the 
institutional review board was obtained prior to imaging. The subject's head position was 
fixed during acquisition using a bite-bar and forehead rest. Pupil dilation was not needed for 
imaging. Scanning areas of the retina were 3x3 mm 2 and 1.5x1.5 mm 2 , as shown in Figs. 2 (a) 
and (b), respectively. Multiple B-scans at the same position, referred to as a BM-scan, were 
used to acquire the phase-change measurements required for the phase -variance contrast. In 
the 3x3 mm 2 scanning pattern, spacing between successive A-scans was 6 um and 15 um 
between BM-scans. Five B-scans were acquired within each BM-scan, resulting in four phase 
changes for the variance calculation at each location. For the 1.5x1.5 mm 2 scanning scheme, 4 
urn spacing between both consecutive A-scans and BM-scans was used. Due to the spatial 
oversampling between consecutive BM-scans, phase changes from within two adjacent BM- 
scans were used for the variance calculations. Three B-scans were acquired within each BM- 
scan in this case, producing two phase change measurements for each BM-scan. For this scan, 
phase-variance calculations utilize four phase changes from two adjacent BM-scans, allowing 
for identical statistics to the 3x3 mm 2 scan pattern. The A-line exposure time was 6.3 us and 
the spectral data were saved in a tiff format for post-processing in Lab VIEW software using 
standard Fd-OCT procedures including spectral shaping, zero-padding, and dispersion 
compensation [31-33]. All retinal images shown in this manuscript were acquired in vivo at 
125,000 axial scans (A-scan) per second. 




Fig. 2. Fundus color photograph of a normal human subject with OCT scanning locations, (a) 
Scanning 3x3 mm 2 of a parafoveal region, (b) Scanning 1 .5x1 .5 mm 2 of a perifoveal region. 

As shown in Figs. 3 (a) and (b), the Fourier transform of spectral fringes (after remapping 
to k-space) from one B-scan acquired by the CMOS detector produces a complex OCT signal 
with intensity and phase information, respectively. Each B-scan consists of 500 A-scans 
acquired over 3 mm lateral scanning range. The dynamic range of the intensity image is 23 dB 
and the phase values vary from -jt to + Jt. Note that one can discriminate retinal vessels, round 
shaped bright structures, in the ganglion cell layer from other retinal structures on the OCT B- 
scan in Fig. 3 (a). However, since the nerve fiber layer is highly scattering as well, a simple 
maximum projection view of the retinal layers doesn't provide sufficient contrast values for 
the retinal vessels. In the Fd-OCT, any particle movement will affect phase values of the 
acquired signal. The basic concept of motion detection relies on calculating phase differences 
of the same locations between consecutive B-scans by 

AO,, (je, z, t) = <D f+1 (x, z,t + T)- O, {x, zj) ( 1 ) 

which can be interpreted as flow induced phase changes in the sample. Here, <S>.(x,z,t) is the 
phase value in the B-scans at lateral location x, depth position z, time t, and T is the time 
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interval between two consecutive B-scans; i is an index of the ;-th B-scan within one BM- 
scan. Then, the following variance calculation identifies mobility without quantitative flow 
information. 

2 tr -i 2 
PhaseVariance = V [AO, (x, z, t) V A€> (x, z, t)f (2) 

where represents the number of B-scans used in the calculation. The preprocessed OCT data 
sets as seen in Figs. 3 (a) and (b) are then postprocessed by software developed at the 
Biological Imaging Center at the California Institute of Technology to extract phase -variance 
signals as shown on Fig. 3 (d). The cross-correlation function between consecutive intensity 
images aligns each B-scan from a BM-scan in an axial direction. Average intensity images 
and phase differences are calculated from sets of B-scans acquired at the same location. The 
image presented in Fig. 3 (c) results from averaging five B-scan intensity profiles, where one 
of the five B-scans is shown in Fig. 3 (a). Phase unwrapping is applied after obtaining the 
phase difference values. Additionally, bulk motion removal and histogram-based thresholding 
processing are implemented to remove phase shifts caused by eye motion. The phase -variance 
based motion contrast procedure generates B-scans of retinal motion information illustrated in 
Fig. 3 (d). The phase-variance data explained above is processed by an algorithm developed in 
Lab VIEW programming software. The differences between results presented by Fingler et al. 
and this paper are primarily limited to changes in the scan patterns, which optimize the retinal 
scanning speed, and the implementation of a high speed CMOS line camera. Further details of 
the post-possessing algorithm and details are explained in the literature [12]. 




Fig. 3. Fourier domain optical coherence tomography of the in vivo human retina, (a) Intensity 
image of a single B-scan. (b) Phase image of the same single B-scan as (a), (c) Average 
intensity image of five B-scans. (d) Phase-variance processed image using phase data from the 
same five B-scans phase data as (c). The scanning size of the lateral direction is 3 mm and the 
sampling density of the lateral direction is 6 um. Scale bar: 500 um. 

To visualize two- and three-dimensional vasculature networks in retinal layers, the three- 
dimensional phase-variance data sets are segmented from the nerve fiber layer to the outer 
nuclear layer. The en face projection view of the segmented data produces two-dimensional 
vascular perfusion maps which are similar to the ones obtained with fundus FA. The FA 
image (1280x1024 pixels) and the color fundus image shown in this paper were acquired with 
a Topcon (TRC-50IX) fundus camera. Moreover, pseudo-color (red-green-blue (RGB)) 
coding of the processed volumetric data sets was used to better visualize three-dimensional 
vessel networks and its projection view creates two-dimensional color vasculature imaging 
linearly scaled by axial location of the vessels. 
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3. Results and discussion 

3.1. Retinal vascular network 

The retinal vascular flow image shown in Fig. 4 (b) was reconstructed using a projection view 
of the retinal contributions of the phase -variance Fd-OCT data. The subject was imaged with 
FA and pvOCT on his left eye. This Fd-OCT volume was acquired over a 3x3 mm 2 retinal 
area captured in approximately 3.6 seconds. Figure 4 (a) shows the fundus FA image over the 
same region for comparison. Note that the phase-variance OCT image is similar to the fundus 
FA image, showing blood vessel networks of the retina. Arteries, vein branches, and relatively 
larger vessels are equally mapped by the two imaging techniques. pvOCT shows some of the 
micro-capillaries but sampling density in this scanning scheme is not sufficient to visualize all 
of them. High speed (125 kHz) imaging of the system allows shorter acquisition time, 
resulting in around 4.5 ms acquisition for each B-scan. This reduces imaging artifacts caused 
by involuntary eye motion. The last B-scan of this volumetric scan data set was used in Fig. 3. 




Fig. 4. In vivo human retinal vasculature images with a normal subject. The size of imaging is 
3x3 mm 2 , (a) Fluorescein angiography, (b) Projection image of phase-variance OCT retinal 
layers. The numbers of A-scans and B-scans are 500 and 200, respectively. The imaging 
acquisition time is 3.6 seconds. Scale bar: 500um. 



In order to reliably visualize micro-circulation, we used a scanning pattern with more 
densely sampled A-scans. In this pattern, sampling spacing is 4 urn between A-scans; BM- 
scans are acquired over a 1.5x1.5 mm 2 area of the human retina. Examples obtained with this 
acquisition scheme, including pvOCT reconstructed B-scans, are demonstrated in Fig. 5. 
Since the acquisition time of each B-scan is approximately 3.5 ms, we expect there is no 
significant eye motion between these scans. An average intensity image composed of three 
successive B-scans from a single BM-scan is shown in Fig. 5 (a). Moreover, Fig. 5 (b) shows 
the image post-processed with the phase-variance method from these data. The combined 
image with averaging intensity and phase variance, red color-coded is shown in Fig. 5 (c). The 
composite image validates that pvOCT identifies most of perfusion structures compared with 
the average intensity imaging where vessels can be classified by relatively bright structures. 
The video (Media 1) of the composite images is attached. Notice that there is a shadow 
artifact in the current phase-variance method that produces residual phase -variance values 
below vessels. This results in non-zero phase-variance values in retinal layers where non- 
perfusion should be observed including in the retinal pigmented epithelium. 

Images illustrated in Fig. 6 are acquired over 1.5x1.5 mm" of the human perifoveal region. 
We compare maximum projection view of pvOCT in Fig. 6 (b) to the fundus FA in Fig. 6 (a), 
acquired at the same location. Figure 6 (b) shows dense micro-capillary networks, not 
distinctly recognizable in the zoomed fundus FA image. The size of the smallest 
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Fig. 5. Human retinal images of phase- variance OCT data processing over 1.5 mm. (a) Average 
intensity imaging of three B-scans from a single BM-scan. (b) Phase-variance image, (c) 
Combined image of average intensity (a) and red color-coded phase-variance imaging (b). 
Media 1 shows the composite image over the 1.5x1.5 mm scanning area. The imaging 
acquisition time of a single BM-scan containing three B-scans is approximately 10 ms. Scale 
bar: 250 pm. 




(a) (l>) (c) 

Fig. 6. In vivo human retinal micro-capillary network. Image size is 1.5x1.5 mm 2 , (a) 
Fluorescein angiography, (b) Projection view of retinal contributions to pvOCT. The numbers 
of A-scans and B-scans are 375 and 360, respectively, (c) RGB depth color-coded projection 
view of retinal pvOCT data. Upper right shows depth scale bar with red color denoting a top 
layer, green showing capillaries as an intermediate vascular bed, and finally blue represents 
micro-capillaries as a deeper vascular plexus layer. The imaging acquisition time of pvOCT is 
3.5 seconds. Scale bar: 250 um. Media 2 shows C-scan fly through of the three-dimensional 
pvOCT data set. 

resolvable capillary in this acquisition scheme is 10-12 urn. Since pvOCT includes depth 
information of capillary locations, Fig. 6 (c) shows a depth color-coded projection view. The 
pseudo-color (RGB) map codes the linearly-scaled axial location of the vessels. The 
superficial capillary network, here mostly within the nerve fiber layer or the ganglion cell 
layer, is coded as red. In addition, the intermediate capillary plexus, between the inner 
plexiform layer and the inner nuclear layer, is denoted by green. Finally, the deeper vascular 
plexus located in between the outer plexiform layer and the outer nuclear layer is colored as 
blue. The attached gray-scaled video (Media 2) demonstrates depth information of the 
vasculature using the pvOCT data set from the nerve fiber layer to the outer nuclear layer. The 
last B-scan of this volumetric data set was used in Fig. 5. 
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Figure 7 shows three-dimensional visualization of depth color-coded image of pvOCT 
from Fig. 6 (c). This image shows volumetric information of the vessel network at the 
perifoveal region in the human retina. Note again that the red-colored, larger diameter 
capillaries have green colored shadows underneath. This is because the pvOCT method 
creates the shadow artifact (non-zero phase-variance values) below relatively large vessels. 




Fig. 7. Three-dimensional depth color-coded image of pvOCT, using the same data presented 

in Fig. 6 (c). 




Fig. 8. Large field of view stitched pvOCT imaging overlaid on a fundus FA. (a) Fluorescein 
angiography, (b) Depth color-coded imaging with the ten volumes. Total image acquisition 
time for ten volumetric images is approximately 35 seconds. 



To achieve a clinically applicable imaging method, the capability to image a large field of 
view is critical. However, the scanning patterns and the acquisition speed presented in this 
paper do not allow for larger than the 3x3 mm 2 scanning area. This was because of the need to 
maintain high-resolution with acquisition time in a single volume of a few seconds. To 
overcome this limitation, we acquired several volumes centered at different retinal locations. 
Figure 8 shows an example of the large field of view image overlaid on the fundus FA. This 
image was made from ten 3x3 mm 2 pvOCT volumes acquired by shifting the position of LED 
fixation points. After pvOCT processing, each projection image was manually aligned with 
linear translation to create a large field of view mosaic. The depth color-coded imaging in Fig. 
8 (b) approximately segmented vessels from below the inner limiting membrane to the inner 
nuclear membrane as a green color as well as from below the inner nuclear membrane to the 
outer limiting membrane as a red color. Moreover, vasculature located in the intermediate 
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layer or laid across two layers segmented above has a yellowish color. The total time of ten 
images for acquisition only is 35 seconds. 

4. Conclusions 

We demonstrate feasibility of a high-speed phase-variance OCT system for two-dimensional 
visualization of the capillary network over foveal and perifoveal regions, as well as three- 
dimensional representation of the vasculature in the human retina. En face visualization of the 
human foveal area with pvOCT produces details of capillary networks comparable to fundus 
FA imaging without using any contrast agents. In addition, densely scanning a small field of 
view, 1.5x1.5 mm 2 , generates the microcapillary perfusion map and color-coded depth 
information of the microvasculature. pvOCT with rapid image acquisition at an A-scan speed 
of 125 kHz reduces acquisition time, hence decreasing eye motion artifacts and producing 
larger field of view imaging. To demonstrate the feasibility of this technique for clinical 
applications, we showed the stitched vascular map of the retina acquired from multiple 
volumes. Our data suggest that pvOCT has potential in the early diagnosis of retinal vascular 
diseases including diabetic retinopathy and vascular -related macular degeneration. As further 
improvements, shadow artifacts will be minimized and high-resolution microcapillary 
networks of larger scanning areas will be visualized. 
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